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Sequence of phase transitions induced in an array of Josephson 
junctions by their crossover to 7r-state 
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PACS 74. 81. Fa - Josephson junction arrays and wire networks 

PACS 64 . 60 . De - Statistical mechanics of model systems (Ising model, Potts model, field-theory 
models, Monte Carlo techniques, etc) 

Abstract. - We show that the transition of Josephson junctions between the conventional and 
7T states caused by the decrease in temperature induces in a regular two-dimensional array of 
such junctions not just a single phase transition between two phases with different ordering but a 
sequence of two, three or four phase transitions. The corresponding phase diagrams are constructed 
for the cases of bipartite (square or honeycomb) and triangular lattices. 



Introduction. — For several decades arrays of weakly 
coupled superconducting islands have been the subject of 
active experimental investigations pQ for many reasons, 
in particular as a simple model system which allows one 
to study the interplay between fluctuations, frustration, 
disorder and other factors in a more controlled situation 
than in bulk superconductors. However, these studies have 
been restricted to arrays of conventional junctions whose 
energy is minimal when the phases of two superconductors 
are equal to each other. 

The first experimental realization of an old theoreti- 
cal idea [2l|3] about fabrication of so-called 7r-junction 
whose energy is minimal when the phase difference on 
the junction is equal to ir was achieved only during last 
decade by Ryazanov et al. [3] who studied superconductor- 
ferromagnet-superconductor (SFS) Josephson junctions 
and observed a transition from the conventional state to 
the 7r-state taking place with the decrease in temperature 
[3] . The experimental investigation of small arrays of SFS 
junctions started almost simultaneously [5J, but insofar 
has been restricted to very modest sizes [5]. 

Since the fabrication of more sizable arrays of SFS 
junctions is definitely a matter of the nearest future, the 
present letter addresses the question what happens with 
a superconducting array of Josephson junctions when the 
decrease in temperature induces a crossover of the junc- 
tions to the 7r-state. Although one could expect (from the 
evident change of the ground state structure) that this 
induces a single first-order transition between two phases 
with different ordering, our analysis reveals that this is 



never the case and in reality an array experiences in the 
crossover region a sequence of two, three or even four phase 
transitions each of which is related with partial or com- 
plete destruction (or restoration) of ordering. The struc- 
tures of phase diagrams and the natures of these transi- 
tions are established both for bipartite lattices (square and 
honeycomb) and for a triangular one. 

Model. — An array of identical SFS junctions can be 
described by the Hamiltonian 



CM') 



(i) 



where tpj is the phase of the superconducting order pa- 
rameter on jth superconducting island, the summation is 
performed over all pairs of neighboring islands connected 
by a junction and V(9) is a periodic even function of 9 
which can have minima both at 9 = and 9 = ir. When 
the contacts forming a junction have low transparency, one 
can keep in the Fourier expansion of 



V{9) = -]T J p cos(p9) 



P =i 

only the first term because a typical value of J p is strongly 
suppressed with the increase of p [3] . 

However, in a SFS junction of an appropriate width the 
decrease in temperature T may force the value of J\ to pass 
through zero and change sign [7] . This leads to the transi- 
tion of the junction from the conventional state [in which 
the deepest minimum of V{9) is at 9 = 0] to the 7r-state (in 
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which the deepest minimum is at 9 = it). Naturally, in the 
vicinity of To, the temperature at which J\(T) = 0, one 
has to keep also the next term in the Fourier expansion of 
V(0), 

V{9) = - Jx cos 6 - J 2 cos(2<9) . (2) 

In the simplest situation the decrease of T leads to the 
change of Ji(T) from positive to negative, while J 2 {T) re- 
mains positive. Our aim consists in analyzing what phase 
transition (or what sequence of phase transitions) takes 
place in a regular array of identical SFS junctions when 
they experience such a transition to the 7r-state (also 
known as — 7r crossover) . 

Bipartite lattice. — First one has to understand 
what would take place with the change of the sign of J\ 
in the absence of thermal fluctuations. For Ji j2 > the 
minimum of the Hamiltonian (JXJ) with interaction ([2]) on 
any lattice is achieved when all variables <p$ (defined mod- 
ulo 27r) are equal to each other, tp- } = Therefore, the 
ground state is characterized by U(l) degeneracy related 
to the simultaneous rotation of all phases. 

The form of the ground state at J\ < depends on the 
structure of the lattice. We start by considering the case 
of a bipartite lattice (square or honeycomb) and after that 
will discuss the more complex case of a triangular lattice. 
For any bipartite lattice the problem with Ji < can be 
mapped onto the problem with J\ > just by rotating half 
of the variables tpj by n. In particular, this immediately 
defines the form of the ground state at J\ < 0, which has 
the same U(l) degeneracy as at J\ > but a different 
(two-sublattice) structure. 

When J 1 = 0, the energies of these two states are equal 
to each other, as well to the energy of any state in which 
all variables ip$ are equal either to $ or to $+7r. Therefore, 
in the absence of thermal fluctuations the system would 
experience at J\ = a single phase transition between the 
phases with different ordering. Note that this property is 
not the consequence of keeping only two terms in Eq. 
- for a more complex form of V(0) the transition will be 
shifted from the point where J\ — to the point where 
the two minima of V(9) have equal depths. However, it 
turns out that in the presence of thermal fluctuations the 
single-transition scenario does not survive. 

The finite temperature phase diagram of the XY model 
with a modified Berezinskii- Villain interaction whose main 
features are analogous to those of Eq. ^ with Ji , J 2 > 
has been constructed in Refs. [5]. In terms of the SFS array 
problem with interaction (J2J the main conclusions of these 
works (confirmed in numerical simulations of Ref. [9]) can 
be reformulated and generalized as follows. 

When both J\ and J 2 are positive and much larger than 
T, the system is in the phase with an algebraic decay of 
the correlation function 

Ci(ji -j 2 ) = (expi^j! -<Pj 2 )}. (3) 

For brevity we shall call this phase ferromagnetic, al- 
though more accurately it should be called a phase with al- 



gebraically decaying ferromagnetic correlations. But since 
in two-dimensional systems with a continuous order pa- 
rameter the real long-range order is impossible |I0j and 
an algebraic decay of correlations |Ilj is as much as one 
can get, the application of such a shorthand is rather nat- 
ural. In terms of SFS array this phase is superconducting 
and is characterized by a finite superfluid density. 

The decrease of J\ down to J\ ~ T induces a phase 
transition of the Ising type related to the proliferation 
of solitons (a soliton is a linear topological excitation on 
crossing which the phase jumps by 7r). The existence of 
such a transition is especially evident for J 2 = 00 when 
the model defined by Eqs. ((TJ) and (0 is reduced to the 
Ising model with coupling constant Ji, however it exists 
(and has the same nature) also when J 2 is less than infi- 
nite. The proliferation of solitons leads to the replacement 
of the algebraic decay of the correlation function C\ (r) by 
an exponential one. On the other hand, on both sides of 
the transition the superfluid density remains finite, which 
for J 2 < 00 manifests itself in the algebraic decay of the 
correlation function 

C 2 (ji -J2) = (exp2i(cp h - cp h )) . (4) 

It is clear that in the phase with such a behavior of 
C\ (r) and C 2 (r) the role of the order parameter is played 
by exp(2ii^>j) and therefore formally it can be called ne- 
matic. Analogous nematic phase (induced by the prolifer- 
ation of solitons) is expected to exist in thin films of su- 
perfluid 3 He [12 . In the nematic phase of a SFS array, the 
superconducting current can be associated with the mo- 
tion of pairs of Cooper pairs and therefore this phase can 
be identified by studying the periodicity of the persistent 
current in the array with annular geometry penetrated by 
a magnetic flux (the period has to be equal to half of the 
superconducting flux quantum). 

The relevant topological excitations in the nematic 
phase are halfvortices, that is the vortices with topologi- 
cal charges ±1/2 which are the end points of solitons. The 
interaction of these objects is logarithmic and keeps them 
bound in pairs, which allows one to treat solitons as closed 
lines playing the role of domain walls in the Ising model. 
With decrease in J 2 the strength of the logarithmic inter- 
action of half- vortices goes down and at J 2 ~ T it becomes 
too weak to keep them bound in pairs. The phase tran- 
sition related to the dissociation of bound pairs of half- 
vortices is of the Berezinskii-Kosterlitz-Thouless (BKT) 
type. It differs from the standard BKT transition by the 
value of the superfluid density jump, which is larger by the 
factor of 4. In the disordered phase the superfluid density 
vanishes and correlation function C 2 (r) also decays expo- 
nentially 

For J 2 <T the disordered phase is separated from the 
ferromagnetic phase existing at large enough ratio J\ /T by 
the standard BKT transition related with the dissociation 
of pairs of integer vortices (exactly like at J 2 =0). With 
the decrease in the ratio T/J 2 one encounters a tricriti- 
cal point, where this BKT transition is transformed into a 
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Fig. 1: (Color on-line) Schematic structures of phase diagrams 
of SFS arrays with (a) bipartite lattice and (b) triangular lat- 
tice. Ferromagnetic (F), nematic (N), antiferromagnetic (AF) 
and disordered (D) phases are separated from each other either 
by continuous (continuous bold lines) or first-order (dashed 
bold lines) phase transitions. The double line separating AF 
and D phases in (b) stands for the sequence of BKT and Ising 
transitions with very small separation. Curved arrows going 
from right to left show different paths of the evolution of an 
array with the decrease in temperature. 

first-order one (with larger than universal jump of the su- 
pcrfluid density). A change in the nature of the transition 
can be associated with switching on of a different mecha- 
nism for the destruction of the ferromagnetic ordering. On 
the other side of the tricritical point the disordering is trig- 
gered not by the integer vortex pair unbinding but by the 
proliferation of solitons taking place when the logarithmic 
interaction of half- vortices is too weak to keep them bound 
in pairs. This induces the simultaneous unbinding of in- 
teger vortices which takes place not because their direct 
logarithmic interaction is insufficiently strong but because 
it is screened by the presence of free half- vortices. 

The schematic structure of the phase diagram contain- 
ing ferromagnetic (F), nematic (N) and disordered (D) 



phases is shown in Fig.QJa) in coordinates J\/T and T / J?- 
Although the above analysis refers only to the right half 
of this figure (with J\/T > 0), in the case of a SFS array 
with a bipartite lattice it is clear from the symmetry of 
the problem that at negative values of J\/T the phase 
diagram has exactly the same form as at positive, the 
only difference being that the phase with the ferromag- 
netic algebraic correlations is replaced by the phase with 
the antiferromagnetic algebraic correlations (which have 
the two-sublattice structure). 

The evolution of a SFS array with the decrease in tem- 
perature is shown in Fig. HJa) by curved arrows going from 
right to left. From the structure of the phase diagram it 
is clear that when thermal fluctuations are taken into ac- 
count the direct phase transition between the ferromag- 
netic and antiferromagnetic phases is no longer possible 
and is replaced by a finite region containing cither one or 
two intermediate phases. 

In particular, for sufficiently low values of To/ ^(Tq) the 
evolution goes along the path F-N-AF, that is, the ferro- 
magnetic and antiferromagnetic phases are separated by 
the nematic phase, both phase transitions being of the 
Ising type. On the other hand, for sufficiently high val- 
ues of To/ ^2 (lb) the ferromagnetic and antiferromagnetic 
phases are separated by the strip of the disordered phase 
and the phase transitions are either of the BKT type or 
of the first order. For intermediate values of T /J2(T ) 
the evolution has to take place along the path F-D-N-AF 
involving three different phase transitions and if in the re- 
gion where \Ji(T)\ is comparable with T or smaller the 
ratio T I JiiT) changes extremely little (by less than few 
percent), the path F-N-D-N-AF involving four phase tran- 
sitions is also possible, although it hardly can be called a 
typical one. 

Triangular lattice. — In the case of a triangular lat- 
tice the structure of the phase diagram at J\ > is basi- 
cally the same as for a bipartite lattice, whereas at J\ < 
the situation is essentially different. The main reason for 
that is that at negative J\ the structure of the ground 
state is different for small and for large values of J\ \ . In 
particular, for — 9J2 < J\ < the minimum of energy is 
achieved when on each triangular plaquette the phase dif- 
ference on two bonds is equal to n and on the third one 
to zero. It is clear that in any configuration satisfying this 
rule the variables <pj can acquire only two values which dif- 
fer by 7r (for example, $ and $ + n), from where it follows 
that in terms of the nematic order parameter exp(2i<^j) 
the system is perfectly ordered. 

After introducing bimodal variables crj = ±1 (below 
they are called pseudospins) such that 

cxp(i( y 9j) = cxp(i$)<7j , (5) 

one finds that the above-mentioned rule is satisfied as soon 
as each triangular plaquette contains both positive and 
negative pseudospins. This means that the set of the al- 
lowed configurations of pseudospins <7j coincides with the 
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set of the ground states of the antiferromagnetic Ising 
model with triangular lattice (the AFMITL model). The 
number of such configurations grows exponentially with 
the size of the system [13]. The exact solution of the 
AFMITL model |13[[T4] at zero temperature is charac- 
terized by an algebraic decay of the correlation functions 
[T5] . in particular, (oj 1 aj 2 ) oc |ji — j 2 1 1 / 2 . These correla- 
tions have the three-sublatticc antiferromagnetic struc- 
ture, that is are positive when the two pscudospins belong 
to the same triangular sublattice and negative otherwise 
[T5] . From the form of Eq. (|5|) it is then clear that at zero 
temperature Ci(ji — J2) coincides with (a^a^) and there- 
fore has a three-sublattice antiferromagnetic structure. 

At Ji < — 9J2 the ground state of ((T|) has exactly the 
same structure as at J2 = 0. In this state each of the three 
sublatticcs is ferromagnetically ordered but the phases in 
the different sublatticcs arc rotated with respect to each 
other by ±2ir/3 [16] . The full set of ground states is char- 
acterized by a combined U(l) x Z2 degeneracy, where U(l) 
corresponds to the simultaneous rotation of all phases and 
Z2 can be associated with antiferromagnetic ordering of 
chiralities of triangular plaquettes. Thus in the absence 
of thermal fluctuations the phase diagram of a SFS array 
with triangular lattice would incorporate three different 
phases, the phases with ferromagnetic and antiferromag- 
netic ordering being separated by a wide strip of the phase 
with perfect nematic ordering and an algebraic decay of 
antiferromagnetic correlations. 

At finite temperatures the perfect antiferromagnetic or- 
dering existing at J\ < — 9J2 is naturally replaced by an 
algebraic decay of Ci(r), however a finite supcrfluid den- 
sity and the genuine long-range order in staggered chirality 
survive. It is known both from numerical simulations |17j 
and analytical considerations [T8] that at J2 = the dis- 
ordering of the system with the increase in temperature 
takes place through the sequence of two phase transitions 
which are situated very close to each other. The first of 
them is related to vortex pairs dissociation and is of the 
BKT type, whereas the second is related with domain wall 
proliferation and is of the Ising type. It follows from the 
analysis of the mutual influence of the topological excita- 
tions of different types [18] that the same scenario can be 
expected to hold also when J 2 > 0. 

The properties of the nematic phase are influenced by a 
small finite temperature more drastically than that of the 
antiferromagnetic phase. It is known both from the exact 
solutions [13l[T4] and from the mapping onto a solid-on- 
solid (SOS) model [H] that at any finite temperature the 
isotropic AFMITL model is in the disordered phase with 
a finite correlation radius (which diverges when T — > 0). 
This immediately allows one to conclude that at T > 
the nematic phase is characterized by an exponential de- 
cay of Ci(r). On the other hand, spin wave fluctuations 
lead to an algebraic decay of 62(1"). These properties are 
in perfect agreement with those of the nematic phase at 
J 1 > 0, which is no surprise since this is just the same 
phase. Exactly like at J\ > 0, at J\ < the nematic phase 



is characterized by a finite superfluid density and its disor- 
dering takes place via BKT phase transition related to the 
dissociation of halfvortex pairs. One more example of an 
XY model in which the phase transition into a disordered 
phase is related to the dissociation of halfvortex pairs is 
the frustrated XY model with dice lattice and one-third 
of flux quantum per plaquctte |20j . 

Since at — 9J2 < Ji < the nematic phase is charac- 
terized by a finite residual entropy So « 0.323 [13], the 
first-order transition line separating it from the antiferro- 
magnetic phase at finite temperatures is shifted to larger 
values of |Ji| (in particular, at low temperatures it takes 
place at Ji m — 9J2 — 2SoT). Together with what we al- 
ready know about the disordering of the antiferromagnetic 
and nematic phases this allows us to draw the schematic 
phase diagram for the case of a triangular lattice shown 
in Fig. DTId). 

Like in Fig. HJa), curved arrows going from right to 
left show the evolution of the system with the decrease in 
temperature. The four arrows present in Fig. [IJb) corre- 
spond (starting from the lowest one) to scenarios F-N-AF, 
F-D-N-AF, F-D-AF, and F-D-C-AF, respectively. Here C 
denotes the phase with long-range order in chirality and 
vanishing superfluid density which separates AF and D 
phases at sufficiently high values of To/J2(Xb). Like for 
a bipartite lattice, the four-transition scenario (involving 
the path F-N-D-N-AF) is also possible if the region where 
J\ (T) j is comparable with T or smaller is sufficiently nar- 
row. 

Conclusion. — In the present letter we have investi- 
gated what happens with a phase- coherent array of SFS 
junctions when the decrease of temperature leads to the 
crossover of the junctions to the 7r-state. The correspond- 
ing phase diagrams have been constructed for the cases 
of a bipartite lattice (square or honeycomb) and of a tri- 
angular lattice. We have shown that the transition from 
the coherent phase existing well above the crossover to the 
coherent phase existing well below the crossover is never 
direct and these two phases are always separated by one or 
more intermediate phase(s). Naturally, the same approach 
can be used to construct the phase diagrams in the vicinity 
of the second crossover (from the 7r-state back to the con- 
ventional state) if it does exist. We hope that our results 
will stimulate more active experimental investigations of 
SFS junction arrays. 

The Hamiltonian (fTJ) can be also used for the descrip- 
tion of a planar magnet with both bilinear and biquadratic 
exchange in the situation when the biquadratic exchange 
is ferromagnetic. For the case of the antiferromagnetic bi- 
quadratic exchange (J2 < 0) such a system with a tri- 
angular lattice has been investigated by Park et al. [21] . 
Its phase diagram also includes a nematic phase, which 
however has a more complex structure of correlations (a 
three-sublattice one) than the nematic phase discussed in 
this work. 

The author is grateful to Ya. V. Fominov for useful dis- 



p-4 



Phase transitions induced in a junction array by — it crossover 



cussions. This work has been supported by the RF Presi- 
dent Grant for Scientific Schools No. 5786.2008.2. 



REFERENCES 

[1] For reviews see Newrock R. S., Lobb C. J., 
Geigenmuller U. and Octavio M., in Solid State 
Physics, edited by Ehrenreich H. and Spaepen F., 
Vol. 54 (Academic Press, San Diego) 2000, p. 263 and 
Martinoli P. and Leemann Ch., J. Low Temp. Phys., 
118 (2000) 699. 

[2] Bulaevskii L. N., Kuzn V. V. and Sobyanin A. A., 
Pis'ma Zh. Eksp. Teor. Fiz., 25 (1977) 314 [JETP Lett, 
25 (1977) 289]. 

[3] For a review of later works see Golubov A. A., Kupri- 
janov M. Yu. and Il'ichev E., Rev. Mod. Phys., 76 
(2004) 411. 

[4] Ryazanov V. V. et al, Phys. Rev. Lett, 86 (2001) 2427. 
[5] Ryazanov V. V. et al, Phys. Rev. B, 65 (2001) 
020501(R). 

[6] Frolov S. M. et al, Nature Phys., 4 (2008) 32. 

[7] Chtchelkatchev N. M., Belzig W., Nazarov Yu. 
V. and Bruder C, Pis'ma Zh. Eksp. Teor. Fiz., 74 
(2001) 357 [JETP Lett, 74 (2001) 323]; Golubov A. A., 
Kuprijanov M. Yu. and Fominov Ya. V., Pis'ma Zh. 
Eksp. Teor. Fiz., 75 (2002) 709 [JETP Lett, 75 (2002) 
588]; Barash Yu. S. and Bobkova I. V., Phys. Rev. B, 
65 (2002) 144502. 

[8] Korshunov S. E., Pis'ma Zh. Eksp. Teor. Fiz., 41 (1985) 
216 [JETP Lett, 41 (1985) 263]; J. Phys. C, 19 (1986) 
4427; Lee D. H. and Grinstein G., Phys. Rev. Lett., 55 
(1985) 541. 



[9] Carpenter D. B. and Chalker J. T., J. Phys.: Gond. 

Matter, 1 (1989) 4907. 
[10] Mermin N. D. and H. Wagner H., Phys. Rev. Lett, 17 

(1966) 1133; Mermin N. D., Phys. Rev., 176 (1968) 250; 

Hohenberg P. C, Phys. Rev., 158 (1967) 383. 
[11] Wegner F., Z. Phys. B, 206 (1967) 465; Berezinskii 

V. L., Zh. Eksp. Teor. Fiz., 59 (1970) 907 [Sov. Phys. - 

JETP, 32 (1971) 493]. 
[12] Korshunov S. E., Zh. Eksp. Teor. Fiz., 89 (1985) 531 

[Sov. Phys. - JETP, 62 (1985) 301]. 
[13] Wannier G. H., Phys. Rev., 79 (1950) 357; Phys. Rev. 

B, 7 (1973) 5017(E). 
[14] Houtappel R. M. F., Physica, 16 (1950) 425; Newell 

G. F., Phys. Rev., 79 (1950) 876; HUSIMI K. and Syozi 
Y., Prog. Theor. Phys., 5 (1950) 177 and 341. 

[15] Stephenson J., J. Math. Phys., 11 (1970) 413. 

[16] Miyashita S. and Shiba J., J. Phys. Soc. Jpn., 53 (1984) 

1145; Lee D. H., Caflisch R. G., Joannopoulos J. 

D. and Wu F. Y., Phys. Rev. B, 29 (1984) 2680. 
[17] Lee S. and Lee K.-C, Phys. Rev. B, 57 (1998) 8472. 
[18] Korshunov S. E., Phys. Rev. Lett, 88 (2002) 167007; 

Usp. Fiz. Nauk, 176 (2006) 233 [Phys. Usp., 49 (2006) 

225]. 

[19] Blote H. W. J. and Hilhorst H. J., J. Phys. A, 15 
(1982) L631; Nienhuis B., Hilhorst H. J. and Blote 

H. W. J, J. Phys. A, 17 (1984) 3559. 

[20] Korshunov S. E., Phys. Rev. Lett, 94 (2005) 087001. 
[21] Park H. J., Onoda S., Nagaosa N. and Han J. H., 
Phys. Rev. Lett, 101 (2008) 167202. 



p-5 



